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Radio frequency spectra of a trapped unitary ^Li gas are reported and analyzed in terms of a 
theoretical approach that includes both final-state and trap effects. The different strength of the 
final-state interaction across the trap is crucial for evidencing two main peaks associated with two 
distinct phases residing in different trap regions. These are the pairing-gap and pseudo-gap phases 
below the critical temperature T^, which evolve into the pseudo-gap and no-gap phases above Tc- 
In this way, a long standing puzzle about the interpretation of rf spectra for "^Li in a trap is solved. 
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Radio frequency (rf) spectroscopy has proved a valu- 
able tool for exploring physical properties of fermionic 
ultra-cold atoms [ll-Q • Interpretation of rf spectra as to 
extract properties of the single-particle excitations, how- 
ever, is quite involved, in particular as it requires one to 
disentangle effects of the trapping potential and of inter- 
actions in the initial and final states of the rf transition. 

After rf spectra were first proposed in Ref.0] to sig- 
nal the presence of a pairing gap for ultra-cold ^Li Fermi 
atoms at low enough temperature, the experimental em- 
phasis has focused in eliminating trap effects via tomo- 
graphic techniques that select the contributions from spe- 
cific trap portions and in reducing the final-state in- 
teraction between the atom excited by the rf pulse and 
the remaining atoms 0. Spectra for ^°K, which are 
not appreciably affected by final-state effects, have been 
recorded for the whole trap in a momentum-resolved way 
Q, mimicking the angular- resolved photoemission spec- 
troscopy widely used in the context of high-Tc super- 
conductivity. From the theoretical side, the spectra of 
Refs.0, 01 were interpreted by accounting for the role 



11 



played by final-state effects in a homogeneous gas [1 
while those of Ref . 0] require one to consider trap but no 
final-state effects. In all these cases the rf spectra show 
a single feature Q with an asymmetric shape. 

Yet, it remains intriguing that the presence of two 
peaks in the original rf spectra of Ref. 3 at temperatures 
about Tc could not be accounted for in a fully satisfac- 
tory way thus far. In Refs.[l3, [3] the two peaks were 
identified as coming from paired and non-paired atoms 
in different trap regions, but final-state effects were not 
included in these works. When applied to a homogeneous 
system, the same theoretical approach predicts two peaks 
in the rf spectra over an extended temperature range [l3| , 
in contrast with the tomographic data of Ref. Q (see in 
particular the supplemental material therein). In addi- 
tion, the work of Ref.[3l yields a single peak for a ho- 
mogeneous system, in agreement with these tomographic 



data, but gives rise to two peaks in a trapped Fermi gas 
only for imbalanced spin populations. We shall show that 
the puzzle posed by the original rf spectra of Ref. Q can 
be solved when final-state and trap effects are taken into 
account simultaneously. 

In this Letter, we present rf spectra of a unitary 
trapped ^Li gas with improved control of the experi- 
mental parameters with respect to Ref. 01 • At the same 
time, we provide a detailed interpretation of these spectra 
on the basis of the theoretical approach of Refs.jlOl. Ill 



which is here extended by a local-density approximation 
(LDA) to account for the trap and by proper inclusion 
of final-state effects also for non-condensed pairs in the 
superfluid phase. 

The main results of our analysis are as follows: 

(i) We are able to confirm the interpretation given at an 
intuitive level already in Ref.01 for the two peaks in the 
rf spectra above (but quite close to) Tc, as associated 
with free atoms in the outer part of the trap and with 
interacting atoms in the pseudo-gap regime in the inner 
part of the trap. [In the present context, a pseudo-gap 
corresponds to a suppression of the spectral weight about 
zero energy detuning.] We have verified that only one 
peak results instead when final-state effects are not taken 
into account; 

(ii) We extend the analysis below T^. and identify the 
two main structures present in the rf spectra as being 
associated with condensed Cooper pairs (a pairing-gap 
phase) and with non-condensed Cooper pairs (a pseudo- 
gap phase). In this case, a minor structure in the calcu- 
lated rf spectra can be associated with free atoms residing 
in the outer part of the trap (a no-gap phase); 

(iii) We verify through a shell analysis, both above and 
below Tc, that the different peaks present in the rf spectra 
originate from different trap regions. 

The experimental starting point of our new rf spec- 
tra is a deeply degenerate, strongly interacting Fermi 
gas of = 4 X 10^ ^Li atoms in an optical trap. 
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Figure 1: Experimental rf spectra of a trapped Li gas 
near Tc (circles) are compared with theoretical calculations 
(DOS+AL) or (DOS) which include or neglect final-state ef- 
fects, (a) Data near unitarity (822G) reproduced from Ref. 0]; 
(b) New data at unitarity (834G). 



The sample is prepared in a balanced spin mixture of 
the two lowest hyperfine sub- levels |1) and |2). The 
basic preparati on p rocedures are described in our pre- 
vious work 0, lid 17|. The mean trap frequency is 
uj = 27r X 196 Hz [18|, leading to a Fermi temperature 
Tp = {ful} /ks) {'iNY^'^ = 1.0 /.tK and a corresponding 
Fermi wavcnumbcr kp = (3800ao)~^, where h is the re- 



duced Planck constant, ks the Boltzmann constant, and 
ao the Bohr radius. We set the magnetic field at the cen- 
ter of the broad 834-G Fano-Feshbach resonance. Our 
thermometry is based on measuring an effective temper- 
ature parameter T and an empirical conversion into the 
real temperature T [B, [H-ilj . 

To record the rf spectra we follow the procedure de- 
scribed in Ref.j^l- By a weak 0.5 s pulse, we drive tran- 
sitions from |2) to the empty state |3) and monitor the 
loss of atoms in |2). After the rf pulse no population can 
be detected in |3), which results from fast three-body de- 
cay in a three-component spin mixture 22, [23| . With 
the scattering lengths — > ±cxd and a/ = — 3310ao [24| 
for the initial channel |1)-|2) and the final channel |1)- 
|3), the corresponding couplings are {kpai)~^ = and 
(kpaf)-'^ = -1.15. 

The main differences of the present measurements to 
the early rf spectra in Ref . 0] are a greatly improved con- 
trol of all experimental parameters , better knowledge of 
the resonance position [2J], and a better thermometry. 

We begin our analysis by reconsidering the data shown 
in the central panel of Fig.l of Ref.[3|, corresponding 
to a magnetic field of 822G near resonance. We report 
those data in Fig[lja), where the energy detuning hio 
with respect to the atomic transition is given in units 
of the Fermi energy Ep = ksTp. The original estimate 
Tp = 3.4 Q has been subject to revision, yielding 
Tp ~ 2.7 /xK, because of an improved characterization of 
the original trap parameters and particle number. The 
temperature at which the data were taken is estimated to 
be 0.3Tf, which is close to T,, 0, while [kpaj)-^ = 0.07 



Figure 2: Shell analysis of rf spectra for a trapped unitary 
Fermi gas at Tc, (a) with and (b) without final-state effects 
when ikpaf)-'^ = -1.15. The regions I, 11, and III are defined 
in the text. 



and (kpttf)^^ = —0.7. This accounts for the differ- 
ences between FigUJa) and Figlljb) which reports the 
new data. 

The double-peak structure in the rf spectrum of 
FiglIJa) was already regarded in Rcf.j^ as an indication 
of the "coexistence of unpaired and paired atoms" . Here, 
we provide the desired quantitative explanation for the 
data of FigllJa), by calculating the rf spectrum above Tc 
extending the theory of Ref. ll[ to account for the con- 
tributions from different shells in the trap through LDA. 
Calculations are presented both in the absence (broken 
line) and presence (full line) of final-state effects. In 
both cases, the theory contains no adjustable parame- 
ters. The first calculation corresponds to a density-of- 
state (DOS) approximation for excitations in the con- 
tinuum which includes pairing in the initial state via a 
self-energy contribution, and the second one to the fur- 
ther inclusion of final-state effects by Aslamazov-Larkin 
(AL) type processes [ll| where the interaction between 
the excited atom and the mate left behind is taken into 
account. 

In this and the following figures, the vertical scale is set 
by normalizing the theoretical spectra to a unit area (sum 
rule) [lo| , while a corresponding normalization cannot be 
exploited for the experimental spectra owing to uncer- 
tainties in the large-w tails. We then make the height of 
the right experimental peak to coincide with the theoret- 
ical prediction. 

The DOS calculation (broken line), which neglects 
final-state effects, cannot reproduce the broad right peak 
present in the experimental data, which is of most in- 
terest being associated with nontrivial pairing effects. 
Within the DOS approximation, pairing effects lead only 
to an asymmetric peak near lo ~ with a slow decay 
oc oj"'^/^ for large w. It is only when final-state effects are 
included in the DOS-I-AL calculation that the presence 
of the second peak at a finite frequency can be accounted 
for, both for its position and width. 

It was shown in Refs.[l3, 111! that for a homogeneous 
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system final-state effects pull the oscillator strength of 
the rf transition toward threshold (w = 0) and make 
the tail to decay faster oc cj~^/^ (the total area of the 
spectrum remaining constant). In a trap, final-state ef- 
fects are stronger where the density is larger, due to the 
larger values of the local final-state coupling {kF{r)af)^^ 
and the smaller values of the local reduced temperature 
T/Tp{r). Final-state effects therefore pull the oscillator 
strength toward threshold, in a more marked way in the 
trap center where the pseudo-gap is larger, and in a less 
marked way when approaching the trap edge where the 
pseudo-gap is progressively reduced. The net effect is 
that the inner and intermediate regions of the trap (I 
and II in the shell analysis of Figl2l^a)) produce partial 
peaks at essentially the same frequency in the rf spectra. 
On the other hand, in the absence of final-state effects 
there is no pulling toward threshold which can compen- 
sate for the spread of pseudo-gap values associated with 
different trap shells, with the net effect of producing only 
a broad shoulder in the total rf spectra (cf. Figl2l(b)). In 
both cases, the peak at zero detuning results from free 
atoms which reside in the outer part of the trap (region 
III in FigslSJa) and[2jb)). In the shell analysis we have 
taken regions I, II, and III to correspond to r/Rp < 0.4, 
0.4 < r/Rp < 0.8, and r/Rp > 0.8, in the order, where 
Rp = [2Ep/{mu]'^)Y/'^ is the Thomas-Fermi radius [H]. 

This interpretation is supported by new data exactly 
at resonance near Tc, which are reported in Figlljb) to- 
gether with the theoretical calculations in the normal 
phase, with and without final-state effects. Note that 
the two theoretical peaks are less pronounced in FiglTJb) 
than in FiglTJa). This is because the final coupling 
{kFaf)~^ = —1.15 of Fig (TJb) is consistently weaker than 
(kpaf)-^ = -0.7 of FigHKa). 

Comparison between theory and experiment is again 
rather remarkable, the only discrepancy appearing in the 
height of the free-atom peak which is somewhat decreased 
in the experimental data (27| . Comparison at higher tem- 
peratures yields further very good agreement between ex- 
periment and theory [28 1. 

We also obtained new experimental data for the su- 
perfluid phase below Tc, as reported in Fig [3] for T = 
O.lTp. Three distinct calculations are reported in this 
case: (a) The DOS approximation extended to the su- 
perfluid phase (broken line); (b) The corresponding inclu- 
sion of final-state effects for the condensed pairs through 
the BCS-RPA approximation of Ref. [lo| (broken-dotted 
line); (c) The further inclusion of final-state effects for 
the non- condensed pairs through an extension of the AL 
contribution of Ref.[l3l to temperatures below Tc (full 
line). At this temperature, the order parameter is non- 
vanishing for r < = 0.63Rp. The DOS approxima- 
tion is used for both r < ro and r > ro, the BCS-RPA 
approximation is included only for r < ro to account 
for final-state effects associated with condensed Cooper 
pairs, while the AL contribution that accounts for final- 
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Figure 3: Experimental rf spectra of a trapped ®Li gas below 
Tc (circles) are compared with theoretical calculations with 
(BCS-RPA, AL) or without (DOS) final-state eflfects. 



state effects associated with non-condensed Cooper pairs 
is included both in the normal phase for r > tq and su- 
perfluid phase for r < ro. 

The agreement between the experimental data below 
Tc and the most complete theoretical calculation (full 
line) appears quite good, in particular for the overall 
shape, width, and wings. The theory predicts the pres- 
ence of a peak at about 0.2Ep and of a shoulder at about 
0.5Ep. This structure is also clearly visible in the exper- 
imental data, with a slight discrepancy in the position 
of the main peak which is shifted to about 0.3Ep. This 
discrepancy could again result from a partial loss of the 
free-atom contribution in the experimental data pTj . As 
a consequence, the theoretical spectrum has more weight 
than the experimental one near zero energy. Note the 
slow decay of the tail of the DOS calculation (iit . 

A shell analysis again helps in understanding the ori- 
gin of the features in FiglS] Accordingly, we report in 
Fig m the shell analysis of the rf spectra below Tc, using 
the same three regions of Figl^l The three panels of FigH] 
correspond to the three theoretical approximations spec- 
ified in the caption of FiglH This analysis shows that: (i) 
The contribution from the inner part of the trap (region 
I) is essentially the same in panels (b) and (c), thus show- 
ing that only a small number of non-condensed pairs can 
be present there; (ii) In the intermediate region II of the 
trap final-state effects are important both for condensed 
and non-condensed Cooper pairs; (iii) The contribution 
from the outer part of the trap (region III) remains the 
same within the three approximations, because in this 
region final-state effects are irrelevant. 

This shell analysis shows that the two main structures 
present in the total spectrum obtained by the most com- 
plete calculation (full line in Fig|4l[c)) correspond to non- 
condensed Cooper pairs residing in the intermediate trap 
region II (left peak) and to condensed Cooper pairs re- 
siding in the inner trap region I (right shoulder). A third 
structure near zero energy is also noticed in the full spec- 
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Figure 4: Shell analysis of rf spectra of a unitary trapped 
Fermi gas below Tc (for (kpaf)'^ = —1.15): (a) Without 
final-state effects; With final-state effects (b) only for con- 
densed pairs and (c) also for non-condensed pairs. As in Fig[2] 
the curves refer to different regions in the trap. 



trum of Figlll^c), which results from free atoms in region 
III of the trap corresponding to a no-gap phase. 

With the support of the shell analyses of Figl2] (above 
Tc) and of Fig|4] (below Tc), we conclude that the trap 
brings about a favorable circumstance: Distinct phases 
are present at the same time in different regions of the 
trap, and final-state effects render them observable as 
distinct features in the rf spectra. 

In conclusion, we have presented an extended compar- 
ison between a new set of experimental rf spectra of a 
trapped unitary ^Li gas, both above and below Tc, with 
improved control of the experimental parameters, and the 
corresponding theoretical calculations that include final- 
state on top of trap effects. These two effects conspire 
constructively with each other, and provide the solution 
to a long standing puzzle about the interpretation of rf 
spectra for ^Li in a trap. They evidence the simultane- 
ous presence of distinct phases, above as well as below 
Tc- The pseudo-gap and no-gap phases above Tc, the 
pairing-gap, pseudo-gap, and no-gap phases below Tc- 
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[27] We speculate that the free-atom peak is more sensitive 
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the peak height in Fig. 1(b) and also lead to the up-shift 
of the peak in Fig. 3. 

[28] See the "supplemental material" . 



Supplemental material 

Comparison at high temperatures 

To test the evolution of the pseudo-gap phase and its 
fate with increasing temperature, we report in FigsISl^a) 
and[SJb) the rf spectra taken at temperatures 0.6Tp and 
1.2Tp, respectively, with the corresponding DOS and 
DOS-I-AL calculations appropriate to the normal phase. 
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Figure 5: Temperature evolution of the rf spectra of a trapped 
unitary ^Li gas. Experimental spectra (circles) for (a) T = 
O.GTp and (b) T = 1.2Tp are compared with theoretical calcu- 
lations which include (DOS+AL) or neglect (DOS) final-state 
effects. In (c) the temperature evolution of the theoretical 
spectra evidences the progressive disappearance of the two- 
peak structure for increasing temperature (here, the spectra 
were rescaled to share the height of the left peak). 



The slow frequency tail of the DOS calculation 



A comment is in order about the comparison which re- 
sults in Fig. 3 of the main text, between the experimental 
data and the DOS calculation (that neglects final-state 
effects). At a first sight it may, in fact, look like that the 
DOS calculation by itself could account for the overall 
features that are present in the experimental data. 
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Figure 6: The dashed line represents the DOS calculation of 
Fig. 3 of the main text, while the circles are the corresponding 
experimental data that have been rescaled vertically as to 
match the theoretical curve near O.SSf. 



However, the comparison reported in Fig|6]shows that, 
if one tries to match the theoretical DOS curve and the 
experimental data near the main structure at O.bEp by 
rescaling the experimental data, one ends up with a tail 
of the DOS calculation (which decays rather slowly like 
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) that considerably overestimates the wing of the 
experimental data. In addition, the DOS calculation 
misses the structure of the main peak at low energy. 

A slow frequency tail resulted already in the theoreti- 
cal approaches of Refs.[Sl,S2] which did not include final- 
state effects. It is only with the subsequent inclusion of 
final-state effects [S3,S4] (which make the tail of the rf 
spectra to decay faster like a;~^/^ and pile up the oscil- 
lator strength toward threshold at the same time) that 
one is able to reproduce not only the main structure at 
0.5Ef but also the tail of the experimental data, as shown 
explicitly in Fig. 3 of the main text. 



The peak at larger frequencies associated with the 
pseudo-gap phase is seen to disappear for increasing tem- 
perature, while the remaining peak at zero detuning 
looses its pronounced asymmetry toward high frequen- 
cies. To evidence this effect more clearly, we report in 
FigEJc) a theoretical DOS-I-AL calculation for five tem- 
peratures in the normal phase between 0.3Tp and O.BTp, 
where the progressive disappearance of the peak at larger 
frequencies for increasing temperature results evident. 
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